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ABSTRACT 



We consider sparticle decays that violate r lepton number, motivated by neutrino oscillation 
data. We work in the context of the constrained minimal supersymmetric extension of the 
Standard Model (CMSSM), in which the different sleptons have identical masses at the GUT 
scale, and neutrino Dirac Yukawa couplings mix them. We find that the branching ratio for 
decay of the heavier neutralino \2 - *■ X + t fj^ is enhanced when the LSP mass m x ~ m? 1 , 
including the region of CMSSM parameter space where coannihilation keeps the relic x density 
within the range preferred by cosmology. Thus xi ~~ * X + t^/-^ decay may provide a physics 
opportunity for observing the violation of r lepton number at the LHC that is complementary 
to r — > /i + 7 decay. Likewise, X2 —> X + e ± /x =F decay is also enhanced in the coannihilation 
region, providing a complement to /x — > e + 7 decay. 
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1 Introduction 



Mixing between different neutrino ffavours has now been amply confirmed by experiments 
on both atmospheric |IJ and solar [jj neutrinos. The distorted zenith-angle distributions 
observed by Super-Kamiokande provided a 'smoking gun' for atmospheric-neutrino oscillations, 
establishing that they are most likely due to near-maximal — v T mixing. More recently, SNO 
has provided two 'smoking guns' for solar-neutrino oscillations, providing direct evidence for 
near-maximal v e — * z/ M T oscillations j§] through its measurements of the charged- and neutral- 
current scattering rates. 

These observations lead one to expect the corresponding charged-lepton numbers to be violated 
at some level. However, the rates for such processes would be unobservably small if neutrino 
masses were generated by the seesaw mechanism [H and there was no lower-energy physics 
beyond the Standard Model. However, the naturalness of the gauge hierarchy, grand unifica- 
tion of the gauge couplings and the relic density of supersymmetric dark matter all suggest 
that supersymmetry should appear at an energy scale ^ 1 TeV. This suggests that processes 
violating the different charged lepton numbers might be observable in low-energy experiments. 
Indeed, charged-lepton-number violating processes could occur at embarrassingly large rates if 
the soft supersymmetry-breaking masses of the squarks and sleptons were not universal. For 
this reason, it is often assumed that these masses are equal at the grand-unification scale, as in 
the constrained minimal supersymmetric extension of the Standard Model (CMSSM). 

Even in this case, renormalization of the soft supersymmetry-breaking slepton masses would 
occur in the minimal supersymmetric version of the seesaw model for neutrino masses, thanks to 
the Dirac Yukawa couplings of the neutrinos |J. These are active in the renormalization-group 
equations at scales between the GUT scale and the heavy singlet-neutrino mass scale, and are 
not expected to be diagonal in the same basis where the light leptons are flavour- diagonal. This 
scenario provides the minimal credible amount of charged-lepton-flavour violation: it could be 
enhanced by GUT interactions and/or non-universal slepton masses at the GUT scale. 

Many signatures for charged-lepton-flavour violation have been considered in this scenario |7], 
H |, [11], [l^], including \i — > and related decays, r — » /i7 and r — > decays. Any or all 



of these may be favoured by the (near-) maximal mixing observed amongst the corresponding 
neutrino species. Other things being equal, one expects these decays to be relatively large when 
the soft supersymmetry-breaking masses m^j and/or mo are relatively small, as has been borne 
out in specific model-dependent studies. Another possibility that has been considered is the 
decay % 2 — > X + e ± /i =F jl3|, |i4)|, where x is the lightest neutralino, assumed here to be the lightest 
supersymmetric particle (LSP), and X2 is the second-lightest neutralino. It has been argued 
that this decay might have a rate observable at the LHC for certain choices of the CMSSM 
parameters. 

In this paper, we consider the alternative decay \2 — > X+ r± A i=F ||15|| . This has certain theoretical 
advantages over the decay X2 —> X + e ± / u T considered previously, since the feedthrough into 
the charged-lepton sector may be enhanced by larger Dirac Yukawa couplings and/or lighter 
singlet-neutrino masses, as compared to the — v e sector, if neutrino masses exhibit the 
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expected hierarchical pattern, and v T — z/ M mixing is also known to be essentially maximal. 
On the other hand, the decay \ 2 ~ y X + T ± A tT has a less distinctive experimental signature 
than X2 — ¥ X + e ± yU =F : . Both decays should be explored at the LHC and a possible linear e + e~ 
collider, and which mode offers better prospects may depend on the neutrino-mass model and 
the experiment. 

We find that the branching ratio for \2 ~~ ► X + r±: A i=F decay is enhanced when m X2 > m fl > m x , 
where f± is the lighter stau slepton. This occurs in a wedge of the (mi/2, Wo) parameter 
plane in the CMSSM that is complementary to that explored by r — ► /ry. The region of 
CMSSM parameter space where this enhancement occurs includes the region where \ — i 
coannihilation suppresses the relic density Q x , keeping it within the range 0.1 < Q x h 2 < 0.3 
preferred by astrophysics and cosmology, even if m^j is comparatively large. The interest of 
this coannihilation region has been accentuated by the latest experimental constraints on the 
CMSSM, such as rrih and b — > 57 decay, which disfavour low values of m\j 2 - We show that 
the branching ratio for \2 — > X + t fi^ decay may be a large fraction of that for the flavour- 
conserving decay X2 X + A i± A iT - An analogous enhancement is expected for the flavour- 
violating decay X2 —> X + e ± /i =F considered by other authors |~3], H , although the absolute 
branching ratio is expected to be smaller. Nevertheless, this decay may provide another way 
of probing lepton flavour violation in the coannihilation region. 



2 Calculational Framework 

We assume the minimal supersymmetric extension of the seesaw mechanism for generating 
neutrino masses, in which there are three heavy singlet-neutrino states iVj, and the leptonic 
sector of the superpotential is: 

W = NftYJijLjHt - Ef(Y e ) i jLjH 1 + ^N C M^ + ^H 2 H X , (1) 

where Y v is the neutrino Dirac Yukawa coupling matrix, M.^ is the Majorana mass matrix 
for the Ni, the Lj and Hj are lepton and Higgs doublets, and the are singlet charged- 
lepton supermultiplets. The superpotential of the effective low-energy theory, obtained after 
the decoupling of heavy neutrinos is 



W eff = LiH 2 (Yj(M D ) V,) LjH 2 — Ei(Y e )ijLjHi. (2) 
\ / ij 

In the basis where the charged leptons and the heavy neutrino mass matrices are diagonal, 

M v = Yj (M D ) ^ Y u v 2 sin 2 (3 (3) 
where the v = 174 GeV and tan/3 = v 2 jv\. 

As mentioned above, we work in the context of the CMSSM, where the soft supersymmetry- 
breaking masses of the charged and neutral sleptons are assumed to be universal at the GUT 
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scale, with a common value m . In the leading-logarithmic approximation, the non-universal 
renormalization of the soft supersymmetry-breaking scalar masses is by an amount 

M)^-^< + Al)(YX)^ g ^. (4) 

We note that, in this approach, non-universality in the soft supersymmetry-breaking left-slepton 
masses is much larger than that in right-slepton masses when the trilinear soft supersymmetry- 
breaking parameter Aq = 0, as we assume here Q. The pattern of charged- lepton- flavour 
violation induced by renormalization depends on the details of 

In plausible mixing textures, the renormalization of the soft supersymmetry-breaking parame- 
ters at low energies can be understood approximately in terms of the dominant non-universality 
in the third-generation left-slepton mass: 

m o LL = diag(mo, m 2 , i x m„), (5) 

where a typical value of the non- universality factor is x ~ 0.9. Correspondingly, we assume 
there is an off-diagonal 7£ — £ll mixing term in the soft mass-squared matrix: 

Am o LL = I 1 - x ) m o — 2 — ' ( 6 ) 

where is the mixing angle between the second and third generation in the charged-lepton 
Yukawa matrix. For the type of non- universalities introduced in (H), this angle can be quite 
large without entering in conflict with the current bounds for r — > fi'y, though in this case large 
mixing in the 2-3 sector must be combined with a small mixing angle between the first and 
second generation, due to the very restrictive bound in the /i — > cy decay [T7]. This mixing 



leads to lepton-flavour violation ~ sin (20), as long as sin(20) is not too large 

We give below numerical results for sample choices of the parameters (x, 0) that may be rep- 
resentative of the possibilities in specific models. We also show how the results vary as (x, 0) 
are varied. 

In the following, we consider mixing between the left-handed r- and yU-flavoured sleptons, but 
f — e mixing might also be present, or even favoured in some models. In such cases, the results 
would be rather similar to those we present, simply with \i replaced by e in the L T - violating 
decay modes studied. 

We consider the following flavour-violating and -conserving \2 decays: 

X2 -> iiij -> Xtftj > X2 -> ViVj -> xWj (7) 
X2 -> XZ -> xti*i 1X2 -> XZ -> XWi (8) 

X2 ^Xh^ xKK (9) 

1 In the case A ^ 0, this parameter would also be renormalized analogously to to? 

2 We have checked that this parametrization is generally a very good approximation for r — ► ^7 decay, as 
well as for \2 X T± ^ ■ 
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The first two decays are the only ones in which flavour violation may be expected, and it would 
of course be unobservable in \2 — > X y v decay. The intermediate sleptons are produced on-shell 
if they are lighter than the \2i while the Z and the h are always on-shell for the range of 
parameters that are of interest to us. Slepton exchanges and h decays may give significantly 
different rates for the various flavor-conserving decays X2 — ► X^t^i ■> suppressing the cases 
£ = /i, e relative to the case £ = r, an effect we see in subsequent plots. 



Our calculations are similar to (|TJj], except that we also include the Yukawa interactions, which 
are relevant for decays into r leptons at large tan/3. Furthermore, we include finite-width 
effects in our calculations of Z° and slepton exchanges. The neutralino and slepton widths, 
which arise mainly from two-body decays, were calculated using the ISA JET package RI8| , and a 
check with calcHEP [l9j found good agreement. For the decays \2 — > X + / i± / i=F ; we found good 



agreement between our code and calcHEP, once we incorporated the VEGAS adaptive Monte 
Carlo programme for the momentum integrals in three-body decays. The results from VEGAS 
differ by several orders of magnitude from those obtained using ISA JET for \2 —> X + / i± A i=F 
decay close to the f resonances. 

For the decay \2 ~ * X + t^t^ the channels mediated by Higgses are important in the areas 
where m X2 — m x < m fl . The widths have been obtained using calcHEP, after adding to the 
package the one- loop QCD corrected Higgs widths from HDECAY 20[| . For flavour- violating 



decays, our calculation agrees with calcHEP, once we modify the MSSM Lagrangian included 
in this package to allow flavour-mixing among fi, r 2 and p,L- 



3 Numerical Results 

The solid (black) lines in Fig. [l| denote the total \2 decay width, as well as the partial widths 
for the flavour- violating and flavour-conserving decays, for the particular cases (a) tan/? = 
10, /i > 0,m 1 /2 = 600 GeV and (b) tan/? = 40, /1 > 0,mi/ 2 = 600 GeV. In both plots, we 
make the representative choices x = 0.9 and </> = 7r/6. In Fig. |l|(a), we see a first edge in 
the flavour- violating width T(x2 — > X + r± A i=F ) & t m o ~ 280 GeV, which is less pronounced in 
T(x2 - ► X + A i± / i=F ) an d almost absent in T{\2 — * X + t ± t =f ). This reflects the dominant role 
of ?2 ~ fi exchange in the flavour-violating case. We also note a second edge when m fl = m X2 
at m ~ 430 GeV, which is visible in all the flavour-violating and flavour-conserving decays to 
X and leptons. The differences between T(x 2 — ► X + t ± t =f ) and T(x 2 X + A i± A iT ) are due, at 
smaller m , to the different masses and couplings of the f 12 and jj>L,R being exchanged, whilst 
the differences at larger m are due to X2 —> X + h decay. 

We see in panel (b) of Fig. [I] features at m = 300 GeV, 420 GeV and 580 GeV, corresponding to 
m X2 = m fl) mf lR and m f2 , respectively. The lowest and highest features show up in T(x 2 — > X + 
r ± /i T ) and r(%2 — ► X+ r±rT ) an d the middle feature in T(% 2 x+/i ± /i T ), as one would expect. 
We note that T(x2 —> X + t ± t =f ) may become relatively large for 300 GeV < m < 580 GeV, 
becoming the dominant X2 decay mode. 

The analogous plot for tan /3 = 10, \i < 0, m = 600 GeV is quite similar to panel (a) of Fig. [I], 
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Figure 1: Comparison of flavour- changing and -conserving X2 decay modes as functions of mo 
for (a) tan/5 = 10,/i > 0,m 1/2 = 600 GeV and (b) tan/3 = 40,/i > 0,mi /2 = 600 GeV. We 
assume for illustration a non-universality factor x = 0.9 and a mixing angle = \ . 



whilst that for tan/3 = 30,/z > 0, m = 600 GeV is intermediate between panels (a) and (b). 
Hence these are representative of the possibilities for flavour-violating \2 decays. 

The ratio of branching ratios R(rfi/fifi) = T(x 2 — > x + r ± /i =F )/r(x 2 — ► X + A i± / i=F ) is shown 
as (red) dashed lines in Fig. |l](a,b). In panel (a), the quantity R(rfi/fifi) also exhibits clearly 
the first edge at m ~ 280 GeV. The second edge at m ~ 430 GeV also appears strongly 
reflecting the facts that flavour violation appears mainly in the left-slepton sector, and that the 
f 2 is mainly f L . We see that, for our choices of x and 0, i?(r/i//i/x) may be of order unity for 
m < 270 GeV, and ~ 10~ 2 for m < 430 GeV. Only at larger m , where the Xi — > X + ^ decay 
becomes kinematically inaccessible, does R(x2 — > X + r± A iT ) drop below 10~ 3 . In panel (b) of 
Fig. [I], we see that R(Tfi/fifi) ~ 0.1 to unity for 350 GeV < m < 580 GeV, dropping below 
10~ 3 only for m > 600 GeV. 

In Fig. H we display contours of the ratio i?(r/x//x/x) of the branching ratios for the flavour- 
violating decay X2 ~ ¥ X + and the flavour- conserving decay % 2 — > x + ^fi T in the x, <p 
plane, for the particular choices tan f3 = 30, m 1( / 2 = 400 GeV and m = 200 GeV of the CMSSM 
parameters. We see that the previous choice x = 0.9,0 = tt/6 is not particularly exceptional. 
To quite a good approximation, i?(r/i//x/i) scales by the square of the factor (1 — x) sin(2</>) 
shown in (^|). This makes it relatively easy to reinterpret our illustrative results in the context 
of any specific flavour texture model that makes definite predictions for x and 0. 

We display in Fig. contours of the branching ratio for the flavour-violating decay r — > ^7 
(thin blue lines) and the flavour-violating ratio R(rfi/fifi) (thick black lines) in the (m^jWo) 
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Figure 2: Contours of the ratio R(ry/ 'pLft) of the branching ratios for the flavour-violating and 
flavour-conserving decays in the x, plane, for tan (3 = 30, fi > 0, m 1 / 2 = 400 GeV and 
m = 200 GeV. 

planes for different choices of tan/3 and the sign of fi. In each case, we have again made the 
representative choices x = 0.9 and = tt/G. 

The contours where R{jfij fifi) = 10 _1 , 10~ 2 , 10~ 3 , 10~ 4 and 10~ 5 are shown as thick black solid, 
dashed, dot-dashed, dot-dot-dashed and dot-dashed-dashed lines. We also display contours of 
BR(t — > yq) = 10 -6 , 10~ 7 , 10~ 8 , 10~ 9 as thin blue solid, dashed, dot-dashed and dot-dot-dashed 
lines. We see that large jx — f (or e — f) mixing is not excluded by the present upper limits on 
BR(t — > fii^e) 1 -/), which are both just above 10~ 6 . We also recall that the X2 is observable at the 
LHC in cascade decays of heavier sparticles EH] for many choices of CMSSM parameters |22 



We see immediately from Fig. |] that the regions where \2 X + r± A i=F ma y De observable 
at the LHC (or a future linear e + e~ collider?), perhaps where R(t/j,/ y/j,) ^ 10~ 2 , are largely 
complementary to those where r — > yy may be observable at the LHC or a B factory, perhaps 
where BR(r -> /i 7 ) > 10~ 8 . 

The darker (green) shaded regions in the bottom right corners of the panels in Fig. are 
excluded because there the LSP is the lighter stau: f\. Such a charged LSP would be in conflict 
with basic astrophysics. The lighter (grey) shaded regions are those in which the cosmological 
relic density of the neutralino LSP x is m the range preferred by cosmology: 0.1 < Q x h 2 < 0.3 as 



calculated using MICR0MEGAS [23], and in agreement with our previous calculations |24L BR |26| 



where h ~ 0.7 is the current Hubble expansion rate in units of 100 km/s/Mpc. In each panel, 
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Figure 3: Contours of R(rfi/ fifi) (thick black lines) and BR(t — > ^7) (thin blue lines) in the 
(mi/2, m ) planes for (a) tan (3 = 10, \i > 0, b) tan (3 = 10, fi < ; (c) tan (3 = 30, /i > 0, (d) 
tan/3 = 40, /x > ; for x = 0.9 and (ft = tc/6. The regions disallowed at low mi/ 2 and m by the 
E821 measurement of at the 2-a level / pg^/ are dark (brown) shaded, the dark (green) shaded 
regions at large mi/2 and low m are excluded because the LSP is the charged f\, the light grey 
shaded regions are those with 0.1 < Q x h 2 < 0.3 that are preferred by cosmology (calculated using 



MICROMEGAS fi25j), and the medium (blue) shaded regions are excluded by b — ► 57 j^/ and the 
dotted line is rrih = 114.1 GeV (calculated using FeynHiggs WH)- 
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there is a region at small mu2 that is disfavoured by laboratory experiments. The regions 
at small m\ii excluded by the b — > sj decay rate are medium (blue) shaded, the regions 
disfavoured by — 2 at small mi/ 2 and m are darker (brown) shaded, and the (dotted) line is 
where = 114.1 GeV as calculated using FeynHiggs p7| . Together, these constraints favour 



the coannihilation strip where ~ m x in all the panels, and the channels at large and 
mo in panel (d) where direct-channel XX ~^ A,H annihilation is relatively rapid. 

In panel (a), for tan/? = 10, \i > 0, the LEP search for the Higgs boson disfavours m 1( / 2 ^ 
360 GeV. In panel (b), for tan/5 = 10,/i < 0, the observed rate for b — > S7 decay excludes 
nii/2 ^ 300 GeV, the LEP search for the Higgs boson disfavours mu2 ~ 430 GeV, and — 2 
excludes a triangle extending up to mi/ 2 ~ 600 GeV. In panel (c) for tan/3 = 30,/i > 0, the 
LEP Higgs limit disfavours m^j ^ 340 GeV, and the other constraints are weaker. A similar 
pattern is repeated in panel (d), for tan/3 = 40, /x > 0. 

In cases (a, b, c), the only region of the (mi/2,mo) plane that survives these constraints is 
the strip parallel to the boundary of the disallowed region, where m x /mf 1 ~ 1.1 — 1.2, and 
coannihilation keeps Q x h 2 within the range allowed by astrophysics and cosmology. This is 
precisely the region where R(rfi/fifi) is maximized, and hence the chances of observing the 
decay may be maximized. We do note, however, that R(rfi/fifi) has a tendency to fall as 7711/2 
increases along this strip, which is apparent in panels (c) and (d). We further note in panel (d) 
that R(t{i/[1{i) ^ 10~ 2 also on the right side of the rapid XX —> A,H annihilation channel, but 
may be significantly lower on the left side of this channel. 



4 Conclusions 



We have demonstrated in this paper that the decay X2 — > X T± ^ T provides an opportunity 
to look for t flavour violation at the LHC that is largely complementary to the search for 
t — > fi'y. Essentially all the above analysis would apply also if the slepton mixing texture 
favours Xi ~ > X r±e=F an d t — > over X2 — ► X T± ^ T an d t — > /ry: it is even possible that both 
Xi X r± / i=F / e=F decays may be observable at the LHC. 

We have phrased this analysis as model-independently as possible. Specific models will predict 
values for the mixing parameters x and 0, and the scaling of our results with these parameters 
is quite simple. We would expect the relevant mixing parameters to be much smaller in the 
case of xi ~^ X/ i±e=F decay, but the corresponding R(/ie/ /x/x) would be enhanced in a similar 
region of the CMSSM parameter space. 

We note that the /i T produced in xi ~^ X r± / i=F decay are likely to have significant transverse 
momentum, and any event in which the X2 is produced is likely to have considerable missing 
transverse energy and jet activity associated with the decays of other sparticles. Therefore, we 
do not expect such events to be suppressed badly at the trigger level at the LHC, though it 
might be more difficult to see X2 ~ > X r±e - T decays. However, a detailed simulation goes beyond 
the scope of this paper. There should be even less problem seeing X2 —* X + T ± A i /e T decays at 
a linear e + e~ collider. We therefore urge more detailed simulations of this decay mode for this 
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machine, as well as for the LHC. 
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